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ABSTRACT 

Two XMM-Newton observations of the fast classical nova V2491 Cyg were carried out in short 
succession on days 39.93 and 49.62 after discovery, during the supersoft source (SSS) phase, yielding 
simultaneous X-ray and UV light curves and high-resolution X-ray spectra. The first X-ray light curve 
is highly variable, showing oscillations with a period of 37.2 minutes after an extended factor of three 
decline lasting ~ 3 hours, while the second X-ray light curve is less variable. The cause of the dip 
is currently unexplained and has most likely the same origin as similar events in the early SSS light 
curves of the novae V4743 Sgr and RS Oph, as it occurred on the same time scale. The oscillations are 
not present during the dip minimum and also not in the second observation. The UV light curves are 
variable but contain no dips and no period. High-resolution X-ray spectra are presented for 4 intervals 
of differing intensity. All spectra are atmospheric continua with deep absorption lines and absorption 
edges. Two interstellar lines of Oi and Ni arc clearly seen at their rest wavelengths, while a large 
number of high-ionization absorption lines are found at blue shifts indicating an expansion velocity of 
3000— 3400 kms^^, which does not change significantly during the epochs of observation. Comparisons 
with the slower nova V4743 Sgr and the symbiotic recurrent nova RS Oph are presented. The SSS 
spectrum of V4 743 Sgr is much softer with broader and more complex photospheric absorption lines. 
The ejecta are extended, allowing us to view a larger range of the radial velocity profile. Meanwhile, 
the absorption lines in RS Oph are as narrow as in V2491 Cyg, but they are less blue shifted. A 
remarkable similarity in the continua of V2491 Cyg and RS Oph is found. The only differences are 
smaller line shifts and additional emission lines in RS Oph that are related to the presence of a dense 
stellar wind from the evolved companion. Three unidentified absorption lines are present in the X-ray 
spectra of all three novae, with projected rest wavelengths 26.05 A, 29.45 A, and 30.0 A. No entirely 
satisfactory spectral model is currently available for the soft X-ray spectra of novae in outburst, and 
careful discussion of assumptions is required. 

Subject headings: novae, cataclysmic variables - stars: individual (V2491 Cyg) - stars: individual (RS 
Oph) - stars: individual (V4743 Sgr) - AAVSO 
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Classical Nova (CN) and Recurrent Nova (RN) out- 
bursts result from nuclear explosions on the surface of a 
white dwarf (WD) that has accreted hydrogen-rich ma- 
terial from a companion star, which in the case of most 
CNe is a low mass main sequence star. Some RNe occur 
in symbiotic systems where the companion is evolved. 
Enough energy is produced to eject the outer envelope. 
The radiative output of the nova first appears primar- 
ily in the optical, but as the density in the ejecta drops 
as a consequence of decreasing mass loss rate from the 
WD surface and continued expansion, the radius of the 
pseudo-photosphere shrinks , and successively hotter lay- 
ers become visible (see, e.g.. lGallagher fc Starrfieldlll978l : 
Hauschildt et al. 19^2)- If nuclear burning continues long 
enough, the peak of the spectral energy distribution even- 
tually shifts into the X-ray regime, and at that time the 
nova emits an X-ray spectrum that resembles those typ- 
ically observed in the c lass of supersoft X-ray sources 
(SSS: Kahabka fc van de n Heuvel 1993 ) . For a review of 
nova evolution see, e.g., TBiidil&_Evani (pool ) . 

In the standard picture of nova evolution, the only time 
at which novae are bright X-ray emitters is during this 
SSS phase, but it is difficult to predict at which time 
this phase starts. In an attempt to catch a nova dur- 
ing the SSS phase, early observations have to be taken 



2 



that bear the risk of no detection. During these cam- 
paigns, several novae have been observed before the SSS 
phase started, yielding faint, hard X-ray spectra (e.g., 
Lloyd ct al. 1992; Kra utter et all 119961 : iMukai fc Ishidal 
2001tlOrio et al.ll2001( ). The early faint, hard X-ray emis- 



sion is believed to arise from shocks in the ejecta or in 
specific cases such as RS Oph and other symbiotic sys- 
tems, in shocks between the ejecta and the dense stellar 
wind of the companion. In very special cases such as 
V458 Vul, early hard emission may also arise from inter- 
actions of the ejecta with surrounding material in a pre- 
existing planetary nebula ([Wesson et al.ll200'8l : lNess et al.l 
[2009a and references therein). 

The uncertainty of the times at which bright SSS emis- 
sion can be seen has been significantly reduced by sys- 
tematic monitoring with the X-ray Telescope (XRT) on 
board the Swift observatory. The first dense Swift mon- 
itoring campaign was carried out after the 2006 out- 
burst of RS Oph ([Bode et al.ll200l lOsborne et aI|[20Tl[: 
iPage et al.1l2008D . and XMM-Newton and Chandra obser- 
vations were sc heduled at times that w ere strategically 
important ('e.g. lNess et"aI]|2007l l2009bV 

V2491 Cyg was disc overed on 2008 April 10.728UT 
by INakano et al.l ([20081 ). The optical brightness was 7.7 
mag at the time of discovery and it reached a peak mag- 
nitude of 7.5 mag one day later. The time scale by 
which the V brightn ess decreased by two magnitudes 
was t2 — 4.6 days ([Tomov et al.l [20088)). as such a 
very fast nova. A rebrightening has been observed af- 
ter JD~ 2454575 (April 26). .Naik et al.. (,2009) dis- 
cussed this secondary peak. In early spectroscopic ob- 
servations, P-Cygni line profiles were found with a ve- 
locity of about —4 000 km s~^ in Ha and H/3 profiles 
(Tomov ct al. 2008b). Near infrared (NIR) p hotometry 
(Naik ct al. 2009) and spectroscopy (jRudv e t al. 2008) 
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obtained in April, 2008 showed no sign of dust formation. 
Unpublished NIR spectroscopy taken in July 2008 also 
revealed no dust formation (Rudy, private communica- 
tion), indicating that V2491 Cyg was not a dust forming 
nova. Rudy ct al. (2008) measured E{B - V) = 0.43, 
which converts to an interstellar neutral hydrogen col- 
umn density of iVn = 2.6 x 10^^ cm~ ^ using the relation 
< A^h/E(B-V)>= 6 ± 2 X 10^1 cm-2 (iDickev fc LockmanI 
[19901: IBohhn et al.l[l978l) . 

The dista nc e wa s estimated 
[Helton et al.' ('2008) and, 

by .Munari ct al. (20l3), both based on a Maxi- 
mum Magnitude versus Rate of Decline rela t ion b y 
Idella Valle fc Livid ((l995l) : [D ownes fc Duerbeckl ([200^ . 
We note, however, that in addition to_the high un- 
certainty of this method (see, e.g., iWarneii I2008D . a 
second, smaller, maximum occurred about three weeks 
after discovery, which reduces the confidence of this 
method. This makes the time zero from which the 
MMRD relation should be used uncertain. 

In X-rays, V2491Cyg was observed with Swift 1.02 
days after discovery, but on ly a marginal detection was 
achieved ([Page et al.l 12010) . which could, however, be 
contaminated by optical loading, i.e., the high number 
of opti cal photons leads to fal se registration of X-ray 
events dib arra fc Kuulkersll2008D . During a second Swift 
observation on day 4.56 after discovery, a count rate of 
0.009 lb 0.002 c ounts per second (cps) was reported by 
iKuuIkers et all ([20081) . In a Suzaku observation taken 9 
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Fig. 1. — Swift X-ray (top) and UV (middle) light curves from 
Page et al. (2010) with times of XMM-Newton observation marked 
by vertical dotted lines. In the bottom panel, the AAVSO light 
curve is added, showing the evolution in the V filter. A small drop 
in X-ray count rate a few days before the secondary peak can be 



days after discovery, superhard X-ray emission extending 
out to 70 k eV was de tected in addition to a hard thermal 
spectrum (jTakei et al. .20091 

This nova is only the second nova to be de- 
tected in X-rays befor e the outburst , after V2487 Oph 
([Hernanz fc Salall2002D . Ilbarra et all ([2009[) reported on 
five Swift pre-outburst observations in which V2491 Cyg 
was serendipitously in the field of view of the XRT with 
clear detections. 

Ba sed on the op t ical sp ectral characteristics of V2491 
Cyg, iTomov et all ([2008aD speculated that it IS a recur- 
rent nova, similar to U Sco and V394 CrA. The similarity 
to V2487 0ph with regards to pre-outburst X-ray emis- 
sion an d the identificati on of V2487 Oph as a recurrent 
nova bv iPagnotta et all (2009) (they discovered a prior 
outburst in 1900) support this speculation. 

The evolution of this no va was followed wit h Swift ob- 
servations every 1-3 days. iPage et all ([2O1O0 give a de- 
tailed description of the evolution of X-ray brightness. 
X-ray spectra, and UV brightness. 

The Swift monitoring observations were used to sched- 
ule tw o XMM-Newton observations as illustrated in 
Fig. [T] ([Ness et aL[|2008al [bl). where, from top to bottom, 
light curves are shown from S'wi/Z/XRT and UVOT (data 
taken from PageNet al. 2010), and V filter magnitudes 
from the database of the American Association of Vari- 
able Star Observers (AAVSO). The secondary peak can 
be seen in the AAVSO light curve around day 15, a few 
days after a small reduction in the X-ray brightness. 

The first XMM-Newton observation was carried out 
five days after the first detection of SSS emission. As 
soon as the indication of a decline was found, a second 
XMM-Newton observation was obtained. A rich dataset 
was acquired, including X-ray and UV light curves, high- 
resolution X-ray spectra in the energy range 0.3-2 keV, 
and low-resolution X-ray spectra in the range 0.1-10 keV. 
In this paper, the soft X-ray spectra and X-ray and UV 



light curves are presented and described. In ^ we de- 
scribe the observations and the data reduction. The hght 
curves are presented and analyzed in ^ and the presen- 
tation of X-ray spectra can be found in 21 These descrip- 
tions need to be complemented by models that will be de- 
veloped and discussed by theoretical working groups, and 
are therefore not dis cuss ed beyond best fits of publicly 
available models in The results from light curve 

and spectral analyses are discussed in fj5] with a focus 
on spectral modeling, comparison with X-ray spectra of 
other novae, and high-amplitude variations. Summary 
and conclusions are given in f|6l 

2. OBSERVATIONS 

Two XMM-Newton observations of V2491 Cyg were 
taken, starting on 2008 May 20.6 (39.93 days after dis- 
covery) and May 30.3 (49.62 days after discovery). The 
XMM-Newton observatory consists of five different in- 
struments behind three mirrors plus an optical monitor 
that all observe simultaneously. The observation details 
are listed in Table [TJ 

For this paper, only the spe ctra from the Reflectio n 
Grating Spectrometers (RGS, I den Herder et al.l I2001D . 
the light curves from the EPIC /pn, and the optical mon- 
itor rOM. lTalaver"al[2009h are used. As will be shown, the 
SSS spectra range from 15 — 38 A (0.33-0.83 keV), a range 
completely covered by the RGS. Since the RGS spectrum 
is sufficiently well exposed for our analysis, the EPIC 
spectra of this spectral range are not needed. At ener- 
gies above 2 keV, optically thin thermal emission is ob- 
servable in the EPIC spectra. These are disc ussed in the 
context of additional Suzaku observations bv lTakei et al.l 
(I20T1I) . 

The RGS consists of two identical grating spectrom- 
eters, RGSl and RGS2, behind different mirrors. The 
dispersed photons are recorded by a strip of eight CCD 
Metal Oxide Semi-conductor (MOS) chips. One of these 
chips has failed in each spectrometer leading to gaps in 
the spectra, fortunately affecting different spectral re- 
gions, and the missing information can be retrieved from 
the respective other spectrometer. A number of spec- 
tral bins contain no information because of bad pixels. 
Again, the redundancy of having two spectrometers al- 
lows almost all gaps to be filled. 

The OM was operated in Science User Defined imag- 
ing plus fast mode, which yields UV light curves. Sev- 
eral shorter exposures were taken with the UVWl (A ^ 
2500 - 3500 A) , UVM2 (A - 2000 - 2600 A) , and UVW2 
(A - 1900 - 2300 A) filters. Some of the OM exposures 
failed because of double bit memory errors (listed with 
exposure time in Table [Ij . These exposures are not 
recoverable. 

In Table [2j the average magnitudes and fluxes in each 
filter are listed for the exposures taken at the times given 
in the first column. The UVW2 magnitudes are consis- 
tent with the Sw ift/VYOT measurements reported by 
iPage et all (|2010l ) for the corresponding times. A graph- 
ical illustration of OM fluxes is shown in Fig. [2 

The EPIC/pn was operated in timing mode. The 
events are all collapsed into a single column, allowing 
fast readout of the chips. While the high time resolution 
is not critical to our analysis, the faster read out prevents 
pile up and the extracted spectra and light curves can be 
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Fig. 2. — OM light curves in flux units. The time axis is plotted 
in units of fractional day after discovery (April 10.728, 2008). The 
time range between days 40.02 and 41.15 are not covered because 
of double bit memory problems leading to unrecoverable loss of 
data. 
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Fig. 3. — X-ray RGSl and OM light curves for the observations 
starting on days 39.93 (top) and 49.62 (bottom). Time bin sizes 
are 50 seconds and 300 seconds for RGSl and OM, respectively. 
The OM count rates were rescaled to flt in the graph. The shaded 
areas mark time intervals over which we have extracted separate 
RGS spectra for studies of spectral changes with the photometric 
variability (see Table [3j. 

used without non-standard corrections. 

We have used standard SAS (Science Analsis Software, 
version 10.0) tools for the reduction of light curves and 
spectra. The new SAS tool xmmextractor was particu- 
larly helpful in determining optimized extraction regions 
and the correction of telemetry losses in the EPIC/pn 
light curve 

3. LIGHT CURVES 

In Fig. [3] we present the X-ray light curves extracted 
from the RGSl in comparison with the simultaneously 

^■^ Many thanks to Aitor Ibarra for support with xmmextractor. 



4 



TABLE 1 
Observation log 



Instrument 


Mode Filter 


Start Time 






End Time 
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ObsID 0552270501 


YYYY-MM-DD@HH:MM:SS 
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OM 
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30 
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TABLE 2 

Ultraviolet magnitudes and fluxes for three OM filters 



UVWl UVM2 UVW2 

A" = 2910 A A" = 2310 A A^^ = 2120 A 

cri cri cri 

Day mag flux^ mag flux^ mag flux^ 



39.93 


12.29 4.39 


± 


0.005 




39.98 


12.30 4.36 


± 


0.013 




40.16 






13.17 2.42 


± 0.01 


40.21 






13.17 2.42 


± 0.01 


40.27 








12.89 3.53 ± 0.021 


40.32 








12.90 3.50 ± 0.021 


49.63 


13.21 1.89 


± 


0.004 




49.68 


13.23 1.85 


± 


0.004 




49.73 






14.06 1.06 


± 0.013 


49.78 






14.03 1.10 


± 0.012 


49.84 








13.79 1.54 ± 0.024 



"effective wavelength 
H0~i''ergcm-2s-iA-i 

taken ultraviolet light curves shown in Fig. The UV 
fluxes are shifted and rescaled to fit in the sample plot. 
The time axis is plotted in units of fractional day after 
discovery (April 10.728, 2008). 

In the first observation, starting 39.93 days after dis- 
covery, the X-ray light curve is highly variable, while the 
UV light curves appear less variable. We have checked 
for variability in the light curves in both bands. For each 
OM exposure, we tested the original (non-scaled) count 
rate, to be constant at the median value and calculated 
a value of reduced x^. AU exposures yield ~ 1) in- 
dicating no significant variability longer than the binsize 
of 50 seconds. The most prominent event in the X-ray 
light curve is an extended dip in count rate, commencing 
on day 40.01, lasting 0.12 days (2.9 hours). After this 



event, the X-ray count rate rises again, yielding a higher 
count rate than before the dip. Unfortunately, two OM 
exposures failed during the most interesting time interval 
between 40.02 and 40.15 days. We are therefore unable 
to study any relations between the dip in X-rays and the 
UV. In addition to the dip, shorter variations of order 15 
min can be seen in the X-ray count rate before day 40.0 
and after day 40.12 which could be oscillations (see §3.11 
below). We fitted a 5th order polynomial to the entire 
RGSl X-ray light curve. Although the polynomial de- 
scribes the long-term trend well, we only found a best-fit 
value of x^ — 10.3 (with 781 degrees of freedom), thus 
additional variability on shorter time scales is present 
at a significant level. Interestingly, during the dip, no 
oscillations seem to be present which we verified by com- 
puting x^ for the time interval during the dip, finding 
X^ < 1 after fitting a 5th-order polynomial to this subset 
as well as for a constant count rate. For the part after 
the dip, we find x^ = 7.2 (340 degrees of freedom) when 
compared to a 5th order polynomial, a bit better, but 
indicating that there is more variability in addition to 
the dip, motivating a detailed analysis of variability in 

The shaded regions indicate three time intervals for 
which we have extracted RGS spectra that are discussed 
in fJD We refer to these episodes as phases a, b, and c 
as indicated inside the shaded areas and the legend (see 
also Table 111). Phase a is the pre-dip spectrum, phase b 
corresponds to the time of low-flux emission, and phase 
c is the post-dip. 

While the OM light curve appears to show some an- 
ticorrelated variations to the X-ray light curves, these 
cannot be considered significant as the OM light curve 
is not variable on a statistically significant level. When 



ignoring the measurement errors, a Spearman Rank test 
yields no correlation or anticorrelation between the X-ray 
and UV light curves. 

During the second observation, the X-ray brightness 
(binned in 50 second bins) deviates only marginally from 
a constant rate, yielding xl ~ 2.9 with 594 degrees of 
freedom. For the OM light curves, binned on 50-second 
grids, the assumption of constant light curves yields xt 
of 1.2 for the UVW2 light curve (which has the lowest 
average count rate), 1.6 for the two UVM2 light curves, 
but 4.2 and 4.6 for the best-exposed UVWl light curves. 
All fits have been done with 42 degrees of freedom. 

The higher degree of variability in X-rays during the 
first observation, especially with the deep dip, is con- 
sistent with a high degree of variations in the early 
Swift /XRT light curve (see Fig. [T|), while the second 
observation was taken at a time when the Swift X-ray 
light curve was also more stable with no dips (jPage et al.l 
|2010( ). Meanwhile, the long-term UV light curve taken 
with the UVOT on board Swift, shows a slow, continu- 
ous decay with no signs of variability (see middle panel 
of Fig. [T|). The XMM-Newton light curves fit into this 
general picture with a high degree of variability in X-rays 
and a lower degree of variability in the UV. 

3.1. Timing analysis 

We have performed period studies using the EPIC / pn, 
EPIC/MOS, and RGS light curves of the first observa- 
tion, all giving consistent results. As discussed above, the 
OM light curves and the second observation are not sig- 
nificantly variable for period studies. In order to remove 
long-term trends, we have detrended the light curves with 
an 8th-o rder polynom ial and applied a Lomb-Scargle al- 
gorithm (Scar gle 1982) and a 2-sinc fitting method devel- 
oped by Dobrotka & Ness (2010). 99. 7% errors are calcu- 
lated from the false alarm probability (jHorne fc BaliunasI 
119861) for the Lomb-Scargle analysis and from S-cr confi- 
dence isocontours for the 2-sine fitting method. Since no 
variability is present during the dip, we have excluded 
the data between 8 and 22 ks of elapsed time. 

The results from the Lomb-Scargle procedure are 
shown in Fig. |4l derived from the detrended light curves 
from EPIC/MOS (black) and EPIC/pn (blue). In the 
top panel, the initial result is shown with two signifi- 
cant frequencies of /i = 0.440 ± 0.033 mHz and /2 = 
0.465 ± 0.035 mHz. The expected location of the first 
harmonic to each detected frequency is marked by a ver- 
tical dotted line, and at most marginal evidence for har- 
monics is present. 

In addition, we have calculated two separate peri- 
odograms for two parts. Part I represents the first 8 ks 
before the dip and part II (22 ks after the dip). The re- 
sults are shown in the bottom two panels of Fig.Hl While 
the first part is too short for any significant detection, the 
second, longer, part clearly yields only a single frequency 
of 0.448 ± 0.036 mHz, corresponding to a period of 37.2 
minutes. 

In order to investigate why the total light curve 
yields two significant frequencies, we have performed 
several tests. W e appli ed the sine fitting method by 
iDobrotka fc NessI (|2010D . and found two frequencies. 
However, a synthetic light curve modulated with a sin- 
gle frequency, sampled and noised as the EPIC/pn data, 
also yields two frequencies. We identify the reason to be 




1.0 

frequency (mHz) 

Fig. 4. — Results from Lomb-Scargle period analysis applied to 
the EPIC/MOS (black) and EPIC/pn (blue) detrended light curves 
while excluding the time interval between 8 and 22 ks of elapsed 
time to avoid contamination by the dip. A false alarm probability 
level of 0.3% is marked by the horizontal line, and all peaks above 
this line can be considered statistically significant at the 99.7% (3- 
cf) level. Two nearby significant frequency values /i and /2 are 
marked in the top panel, together with the expected frequencies 
of their first harmonics. Separate periodograms calculated for the 
first 8 ks of data before the dip (part I) and those after 22 ks (part 
II) have been computed and are shown in the bottom two panels. 

an increase in the amplitude of variations after the dip, 
which mimics a beating cycle when fitted with two sine 
curves. This is illustrated in Fig. [5] where the detrended 
EPIC/pn light curve is shown in grey with the best- fit 
2-sine curve overplotted. One can clearly see that the 
changes in amplitude in the model reproduce the ampli- 
tude change in the data. While the change in amplitude 
could be beating of two signals, we consider it more likely 
that some other process that is related to the dip is re- 
sponsible for the change in amplitude. 

For the discussion of the origin of the oscillations it 
is important to know whether the first harmonic is also 
present. In the bottom panels of Fig. [5l we show the 
first and second parts with two different light curve fits 
based on the frequencies detected in the respective peri- 
odograms (see bottom two panels of Fig. U]). The light 
blue lines are single sine curves which were obtained after 
fitting amplitudes and phase while keeping the frequen- 
cies fixed at 0.381 mHz for the first part and 0.448 mHz 
for the second part. The black curves are fits with fixed 
fundamental frequencies plus the respective first har- 
monic (twice the frequency values), also with variable 
phase shifts and amplitudes. While the detection of a 
first harmonic is in no case statistically significant (as 
can be seen from Fig. [4]), the best fit to the first part 
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Fig. 5. — Top panel: Result from fitting two sine curves to the de- 
trended EPIC/pn light curve. The data are in light grey and model 
in black. The strongest peak in the Lomb Scargle periodograms 
are shown in the last two panels in Fig. |4] and as sine curves over- 
plotted over the first part (middle) and second part (bottom) . The 
blue curves are single frequencies of 0.38 i mHz and 0.448 mHz, and 
the black curves are the sum of the ground frequency and the re- 
spective first harmonic. 

suggests that it may be present but is not significant 
only for tlie reason tliat no more than three cycles are 
covered. We can therefore not fully exclude that a first 
harmonic was present before the dip which then disap- 
peared in the dip, as is clearly not present in the second 
part. 

The frequency at 0.448 ± 0.036 mHz (with a pulse frac- 
tion of 7% after the dip) seems clearly detected, but we 
wish to caution that only seven cycles are covered, bear- 
ing a certai n risk of red noi se. This has been discussed 
in detail bv iVaughanI (|201Cll ) and references therein. We 
tested the c onfidence of our signal using the method 
proposed bv IVaughanI (2005). The 0.448-mHz peak is 
approximately at the 98% confidence limit. More data 
could mitigate this situation, but the entire observation 
cannot be used because the dip may have changed the 
conditions under which a periodic signal is visible. The 
second observation yields no periodic signal. Our con- 
clusions are based on seven apparently periodic cycles in 
the second part of the first observation. 

4. SPECTRA 

The SAS tool rgsproc was used to generate RGS 
events files from the raw data and extract spectra in 0.01- 
A bins from a standardized extraction region around the 
dispersed photons. Time filters can be applied, allow- 
ing for the accumulation of spectra over any given time 
interval. The SAS tool rgsf luxer combines all informa- 
tion from both spectrometers into a single spectrum in 
photon flux units. It does not correct for redistribution 
of monochromatic response into the dispersion channels. 
Thus, the 'fluxed' spectrum is not suitable for quantita- 
tive analyses but it can be used for the purpose of vi- 
sualization and comparison with broad-band models or 
other observations. 

4.1. Description 

In Fig. |6l the fluxed RGS spectra from four different 
epochs are shown. The color code is the same as that 



used in Fig. [31 In Table |31 the respective start times in 
units of day after discovery, exposure time, integrated X- 
ray band fluxes, and best-fit blackbody parameters (see 
below) are listed. The fluxes are absorbed fluxes and 
correspond to logarithmic luminosities of 30.67, 30.11, 
30.82, and 29.98, respectively (in cgs units). The spec- 
tra are all continuous spectra with deep absorption lines 
at, e.g., ~ I9A and ~ 24.5 A. The overall shape of the 
spectra resembles that of an absorbed blackbody, which 
is demonstrated by the overplotted best-fit blackbody 
curves that were obtained by minimization. Spec- 
tral bins around deep absorption lines have been dis- 
carded for better reproduction of the continuum. The 
effects of photoelectric absorption along the line of sight 
were modeled using the Bo uldcr-Tiibingen model TBabs 
developed by I Wil ms et al.l ([2000) with variable oxygen 
abundance. The best-fit parameters are listed in the last 
four columns of Table [3l where the last column lists the 
oxygen abundance in the TBabs model, relative to solar 
(jGrevesse fc Sauvalll998| ). It was allowed to vary in or- 
der to reproduce the depth of the Oi absorption edge 
at ~ 23 A. The corresponding ls-2p absorption line at 
23.5 A is also modeled as part of TBabs. The other ab- 
sorption lines are of photospheric origin, indicating that 
the spectrum is that of an atmosphere. The deepest pho- 
tospheric absorption lines at ~ 19 and ~ 24.5 A have 
been excluded from the fit. 

With the given complex structure of the expanding 
nova ejecta, the blackbody fits are only parameteriza- 
tions. As shall be shown in MA\ obtaining reliable results 
from atmosphere modeling is sufficiently complex to re- 
quire a separate project, and the blackbody temperature 
is here given only as a preliminary characterization of 
the spectral hardness. The statistical uncertainties are 
small, of order 3%, because of the high quality of the 
data, however the values of reduced are 2.4, 1.8, 8.4, 
and 2.0 for days 39.93, 40.03, 40.18, and 49.62, respec- 
tively, and standard error determination is not reliable 
for large values of reduced x^- Blackbody fits without 
discarding any spectral bins yield systematically higher 
temperatures by ~ 0.5 x 10^ K, with much larger values of 
reduced of 6.3, 2.2, 28.8, and 3.5. We emphasize that 
these temperatures cannot be interpreted as the photo- 
spheric temperature because a blackbody model does not 
account for the most fundamental physics in the ejecta. 

For comparison with the observed UV fluxes listed in 
Table [2l we have calculated the projected fluxes from the 
blackbody fits. For days 39.93, 40.03, 40.18, and 49.62, 
we found (0.61, 1.5, and 2.1)xl0-^'^ergcm-2 s'^ A"!, 
(0.06, 0.15, and 0.21) x 10"" ergcm'^ s'^ A"!, (0.43, 
1.1, and 1.5) X 10^1" erg cm-^s-iA-i, and (1.9, 4.9, and 
6.9)xl0-i^crgcm-2s-i A-i for the UVWl, UVM2, and 
UVW2 filters, respectively. These fluxes are significantly 
lower than the measured values. Since other sources of 
UV emission are likely present, these numbers do not rule 
out a blackbody nature of the broad-band spectrum. 

The changes in blackbody temperature are fairly small 
between the three epochs of the first observation, w hich is 
consistent with the findings bv lPage et al.l (|2010| ). How- 
ever, the second observation yields a significantly harder 
spectrum. At the same time, the best-fit value of A^h is 
of the same order as the interstellar value derived from 
E{B—V) measurements, indicating that the previous ob- 
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Fig. 6. — RGS spectra extracted from different epochs. The time intervals over which each spectrum has been extracted are marked 
in Fig. [3] with the corresponding colors indicated in the left legend (see Table [3] for exposure details). The units are flux units, and the 
spectra are combined from RGSl and RGS2 using the SAS tool rgsfluxer. The thin lines are best fit blackbody fits with the blackbody 
temperatures given in the right legend. 



TABLE 3 

Separation of RGS spectra into four time segments with blackbody and atmosphere model parameters 



Start 
Date" 


Color 
Code 


At'' 
ks 


flux'^ 


-'bb 
lO^K 




r>d 
^BB 
10* cm 




-'off 

lO^K 


"h 


cir 

10* cm 


logC^bol)" 
erg s~^ 


Mo 




9 


39.93 


blue (a) 


4.7 


3.57 


5.9 


4.9 


77 


1.6 


9.6 - 10.3 


2.1 - 2.4 


5.7 - 7.3 


38.42 - 38.52 


2.48 - 4.18 


3.4 - 


4.1 


40.03 


pink (b) 


6.1 


0.98 


6.3 


4.1 


22 


2.0 


9.8 - 10.5 


1.6 - 1.8 


2.5 - 3.1 


37.73 - 37.81 


0.46 - 0.72 


4.7 - 


5.4 


40.18 


red (c) 


16.5 


5.09 


6.3 


4.7 


57 


1.2 


9.7 - 10.4 


2.0 - 2.3 


5.8 - 7.8 


38.45 - 38.58 


2.54 - 4.34 


2.4 - 


2.8 


49.62 


purple 


29.8 


0.72 


8.2 


2.7 


4 


2.3 


10.3 - 10.5 


1.5 - 2.0 


1.9 - 2.2 


37.47 - 37.62 


0.27-0.37 


4.4 - 


4.7 



"day after discovery, 2008, Apr. 10.7 
''Net Exposure time 

c]^g-lO gj-g g-l QYgj. range 11-37 A 

''from blackbody fit 
■^from atmosphere model 

/ 1021 cm-2 

s Oxygen abundance in TBabs model, relative to solar bv lGrevesse fc Sauvall 11998) 



servations contained exces s absorption from local neutral 
material. Note also that iPage et all ()2010 ) found even 
lower values of A^'h around day 50 (see their fig. 4) . While 
the increases in spectral hardness could be due to an in- 
creasing effective temperature, which would be consistent 
with a smaller effective blackbody radius, changes in the 
atmospheric structure may also have occurred. Also, an 
increasing degree of ionization of circumbinary neutral 
m aterial can l e ad to spectral hardening as demonstrated 
bv lNess et all (|2007D for RS Oph. The degree of ioniza- 
tion of oxygen influences the depth of the O i absorp- 
tion edge at 23 A and the associated ls-2p line at 23.5 A. 
The parameter A{0), the abundance of neutral oxygen 
in the cold absorber component, yields higher values for 
a lower degree of ionization because more neutral oxygen 
is present. The results given in the last column of Table |3] 
indicate that on day 40.03, when the total flux was lower, 
the degree of ionization was also lower, which can be ex- 
plained by recombination in a photoionized plasma that 
was less irradiated during the times of low- flux emission. 

Some model-independent conclusions can be drawn by 
direct inspection of the details in the spectra. Trends 
of abundances can be estimated from line depths, and 
the local expansion velocity can be determined from line 
shifts on a quantitative level. 



4.2. Details in the High-resolution spectra 

The details of the spectra are shown in Fig. [71 using the 
same color code as in Figs. |3] and [Bl Line labels of promi- 
nent transitions are placed at their rest wavelengths us- 
ing different colors that indicate groups of transitions 
and elements. The groups are: neutral (interstellar) lines 
(black), the H-like Ly series (Is-np) of Nvii, Oviii, and 
Cvi (black, red, and blue, respectively), the He-like Is- 
np series of Nvi, Ovii, and Cv (purple, orange, and 
grey). Neon emission lines (cyan), and absorption fea- 
tures that cannot be identified (green) . Lines that are not 
detected but are of interest are marked in light grey (e.g., 
Cv). Expected blue shifts corresponding to a Doppler 
velocity of —3200 km s~^ are indicated by dotted lines 
for each transition. For the majority of lines, the dotted 
lines coincide with deep absorption troughs. Except for 
the two interstellar lines of Oi (23.5 A) and Ni (31.28 A), 
all lines indicate roughly the same expansion velocity of 
-3200 kms"^ 

In the spectrum representing the low emission phase 
during the first observation (phase b, pink), the Nvii 
and the C VI Lya lines at Aq = 24.78 A and Aq = 33.74A, 
respectively, appear to be blue shifted by a larger amount 
than in the other two spectra from the first observation 
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Fig. 7. — Combined 'fluxed' RGSl and RGS2 spectra of V2491 Cyg extracted from different time intervals. The start times are given in 
the upper left legend, and the colors correspond to those used in Fig. [3] (see Table[3ll. Important absorption lines are marked with small 
vertical solid lines at their rest wavelengths and with dashed lines at the wavelengths corresponding to blue-shifts of ti = — 3200kms~^. 
The line labels are color coded indicating lines of N VII in black, N VI in purple, O IV- VI and O VIII in red, O VII in orange, carbon lines in 
blue, interstellar lines in black, and unidentified lines in green. Labels in grey indicate lines that are not detected. 
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Fig. 8. — Same as Fig. [7] in linear units, focusing on the Ovil 
edge at 16.77 A. 



(phases a and c in blue and red). Meanwhile, the Ly/3 
lines of the same elements and the Nvi lines show the 
same blue shift. 

The H-like Ly series lines of N vii can be identified up 
to ls-5p (LyS). The recombination/ionization continuum 
of N VII is expected at 18.59 A (667.1 eV) but is not seen 
(and is not marked) because it overlaps with the Ovii 
ls-3p line at 18.67 A. For Oviii, the Is-np series can be 
seen up to n = 4 (Ly7). At shorter wavelengths, not 
enough continuum is present to detect higher-n Ly series 
lines of O viii. For C VI, the Is-np series can only be seen 
up to Ly7, while hyd and the recombination/ionization 
continuum at 25.3 A (490 eV) are not present. The He- 
like series lines of Nvi and Ovii are detected up to Is- 
5p, and for Ovii, the recombination/ionization contin- 
uum at 16.77 A (739.3 eV) is present between 16.6-16.8 A. 
This is shown in more detail in Fig. [H The depth of the 
line is strongest in the spectrum from day 40.18. The 
shape of the edge indicates an imbalance in favor of ion- 
ization, thus O VIII is being pumped by photoionization. 
The N VI recombination/ionization continuum at 22.46 A 
(552.1 eV) cannot be detected because it overlaps with a 
number of low-ionization oxygen lines, O vi, O v, and 
O IV, marked in red. While the He-like ls-2p line of C v 
(40.27 A) is outside the spectral range of the RGS, the 
high-order series lines are in the RGS band pass and are 
not detected (grey labels). 

Unidentified absorption features are marked with a 
question mark at a wavelength that would be the rest 
wavelength if that feature was photospheric and thus blue 
shifted by the same amount as the other lines, but they 
could also be interstellar lines. Identification of these 
lines might be possible if they were also seen in other 
nova spectra with a different velocity profile (see i 35.2.2|) . 
Seeing similar features at the same wavelengths would 
be an argument in favor of interstellar lines, while line 
shifts by an amount commensurate with other line shifts 
would be more suggestive of photospheric lines. 

At 12 and 13.5 A, the H-like and He-like hues of Nex 
and Neix are clearly detected in emission. On day 49.62, 
this wavelength region is dominated by continuum emis- 
sion originating from an increased Wien tail of the pho- 
tospheric continuum spectrum. In the logarithmic plot 
it is difficult to see, but the Ne li nes are still signific antly 
present, albeit with a lower flux (jNess "ei~ani2008bn . 

4.3. Line profile modeling 



The deepest absorption line is the Lya line of Oviii 
at 18.8 A (Ao =18.97 A). iNessI (|2010[ ) found that it was 
saturated without going to zero. This is the only line 
for which the bottom of the line profile is flat, indicating 
that it has the largest column density and thus the high- 
est number density. In addition to a high oxygen abun- 
dance, the temperature structure at the time of the ob- 
servation provides particularly advantageous conditions 
for the formation of the O viii ionization state. The lat- 
ter is supported by the presence of the Ovii to Oviii 
ionization edge at 16.8 A (see Fig. [8]), which is an indica- 
tor for photoionization. 

We have u sed a s lightl y modified version of the method 
described bv lNessI ()2010[) to determine hue shifts, widths, 
and optical depths at line center for each of the four 
spectra from days 39.93, 40.03, 40.18, and 49.62. A 
line profile model spectrum is calculated over a narrow 
wavelength ra nge around an absorption line, based on 
equation 1 in iNessI (|2010( ) . For faster computation, we 
perform no correction for interstellar absorption (thus 
T(A) = 1), and for the continuum we assume a lin- 
ear function rather than a blackbody curve. These two 
effects are marginal over the small wavelength region 
around the absorpt ion lines. We have compared the re- 
sults by INessI (120101) with the same spectra and lines but 
with the modified method and get consistent results. 

The resulting parameters of interest are the optical 
depth at line center, Tc, the central (rest) wavelength, 
Ao, and the Gaussian line width, a. The line shifts, 
(A— Ao), and line widths were converted to Doppler veloc- 
ities Wshift = c(A — Ao)/Ao and Uwidth — ca/Xn, where c is 
the speed of light. The instrumental line broadening can 
be estimated as finstr ~ 300 km s""'^ and was accounted 
for in computing Wwidth = ^(ccr/Ao)^ - vf^^^.^. 

In addition to these parameters, line column densi- 
ties, A'x, were obtained by integrat ion over th e opac- 
ity t(A) as defined in equation 2 of INessI ()2010 [). using 
oscillator strengths extracted from the Chianti database 
(pLandi et al.ll2006l ). which is an approximation for strong 
lines. 

All results are summarized in Table HI where the top 
block lists the results for the two brightest spectra taken 
before and after the dip in the first observation (blue and 
red shadings in the top panel of Fig. [3]) , and the bottom 
block represents the two faintest spectra represented by 
pink and purple colors in Fig. [T] Given are name of ion, 
ionization parameter (see last paragraph of this subsec- 
tion), rest wavelength, oscillator strength, line shift and 
width, optical depth at line center, and line column den- 
sity. The given errors were obtained for each parameter 
from a grid of parameter values versus corresponding 
by interpolating the two parameter values on each side of 
the minimum that are larger by 3.53 than the minimum 
X^. This increase in with three free parameters of in- 
terest, yields 1-cr uncertainties. While stepping through 
the grid, all other parameters were iterated, including 
the parameters for the contin uum that are statistically 
'uninteresting' fsee IAvnilll976( ). starting at their best- fit 
values. If a better fit was found during this procedure, 
a new minimization was started with the new parame- 
ters until the fit was stable. With this procedure, we are 
confident to have found global minima in x^- Only in 
one case have we not accepted the best fit, which is illus- 
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Fig. 9. — Line profile fitting to the spectrum of day 40.18, fol- 
lowing a modified method by (Ness 2010, see text) for the two 
blended lines of N VII (Ao = 24.78 A) and Nvi (Aq = 24.90 A). The 
observed spectrum is plotted in histogram style (blue for fitted part 
and black for non-fitted part), the best-fit model with a thick red 
line (solid and dotted for fitted and non-fitted, respectively), and 
the model components for each line with grey dotted lines. While 
the fit in the top panel yields lower x'^, we consider the fit in the 
bottom panel as the one with better physical motivation because 
the parameters are closer to the values from the other spectra (see 
Table lU, and a dip between —2000 and — 1800kms~i is better 
reproduced. 

trated in Fig. ^ where in the top panel the global best 
fit is given while in the bottom panel our preferred fit is 
shown, which is a local minimum in x^- The best fit re- 
quires a broad N VI line and a narrow N vii line, while all 
other observations yield the opposite result. This single 
result does not agree with the other results. We illustrate 
this to make the point that rigorous fitting does not al- 
ways lead to the best results. In Fig. ^ a dip between 
—2000 and — 1800kms~^ can be seen which could well 
be originating from the N vi transition. In the bottom 
panel, this dip is much better reproduced, yielding a ve- 
locity shift for Nvi that is consistent with that of the 
large majority of other absorption lines. Although this 
fit has a higher value of x^, we prefer this fit over the 
one in the top panel, and in Table 01 the best-fit values 
are given. 

Fig. [To] is a graphical illustration of the line shifts 
and line widths (top panel) for different lines at differ- 
ent times. Only strong ls-2p (Lya and Hea) transitions 
were selected from Table H] for this plot for better clar- 
ity. The majority of lines are blue shifted by —3000 to 
— 3400kms~^ and are broadened by up to 1200 km s~^, 
which is significantly broader than the expected line 
width from instrumental broadening. We can thus con- 
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Fig. 10. — Line shifts plotted versus line widths (top) and versus 
ionization parameter yielding the highest column density (bottom) , 
with the values given in Table |4] The symbols indicate different 
lines as indicated in the right legend and four symbol colors corre- 
spond to the four epochs of extracted spectra (Fig. [3l and TablelSj. 

elude that the ejecta are significantly extended, allowing 
us to view a range of expansion velocities and thus a 
range of different plasma layers. 

The Oviii line clearly stands out with consistently 
larger line shifts and widths in all three data sets of the 
first observation. By day 49.6, the line has become nar- 
rower but has the same line shift. 

The Nvii and Cvi lines are blue shifted by a larger 
amount during the dip in the first observation (day 40.03 
in pink). Together with the Oviii line, these are all 
H-like transitions while the He- like Nvi line, which is 
formed at lower temperatures, is not more blue shifted 
during the dip. This could indicate that the high- 
temperature lines probe higher wind velocities which are 
likely found further outside. 

To get a better idea about systematic temperature ef- 
fects, we computed a characteristic ionization parameter 
^ = L/(nr^), with L the bolometric luminosity of the 
ionizing source in the 1 — 1000 Rydberg range, n the 
electron density of the photo-ionized absorbing gas, and 
r the distance from the ionizing source. For a grid of 
ionization parameters, column densities have been com- 
puted for each ion usin g ionizatio n-balance calculations 
done with Cloudy (.Ferland et al. 1998). The spectral 
energy distrib ution (SEP ) input into Cloudy consists of 
IR data from iNaik et all ()2009f), optical data fr om the 
AAVSO database presented by iHachisu k Katol ([2009), 
UV and soft X-ray ffuxes f rom this paper, and hard X- 
ray fluxes from .Page et al.l pOlOi ). Solar abundances by 
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TABLE 4 

Results from line profiles measurements 



Ion log(g) 

(ergs"-*- cm) 



Ao 

A 



■"shift 

kms"-'^ 



""width 

kms^-*- 



IQl® cm" 



'^shift 

kms^-'- 



^'width 

kms^-*- 



TVx 



Day 39.93 

Cvi/3 1.30 
Cvi« 1.30 
Nvi/3 0.80 
NviQ 0.80 
Nviio 1.60 
Ovii/3 1.20 
Oviio 1.20 
Oviii/3 2.70 
O villa 2.70 



28.46 
33.74 
24.90 
28.78 
24.78 
18.63 
21.60 
16.00 
18.97 



0.16 
0.83 
0.14 
0.66 
0.83 
0.15 
0.68 
0.16 
0.83 



-2705ll« 

-soiotii; 

-3157t^« 
-31191^^ 
-3275tiO^ 
-3227li^o 

— 3292"''^^^ 
-3703 ± 107 



< 353 
0"'-191 

631115" 
1146ti^^ 

'^°-200 
97^9+613 
^'0^-493 

2545ti^ 



0.37; 
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Fig. 11. — Comparison between model atmosphere and RGS 
spectrum of V2491 Cyg comm encing on day 40.18. Models are 
taken from'Rauc li et an i l2010bl ). Their table 1 lists different abun- 
dance categories, and the blue solid and red dotted lines are mod- 
els #007 and #003, respectively. Both models are blue shifted by 
z = —0.01, corresponding to tighift = — 3000kms~^. All models 
assume logg = 9. 



ILodders fc Pah^ (1200^ were assumed. A description of 
how t his was done can be found in §6.2 of the SPEX man- 
uaP^ (jKaastra et al.lll996l) . From the resuhing balance 
curves, the ionization parameter yielding the largest ab- 
sorption column density for each ion is selected for plot- 
ting along the vertical axis in the bottom panel of Fig.fTUl 
As a confidence range, the range of ionization parameter 
is chosen that includes all column densities greater than 
10% below the value at the peak. A trend seems sugges- 
tive that lines with higher ionization parameter are more 
blue shifted, but this conclusion depends very much on 
the measurements of the O Viii line. With the given set 
of lines, it is thus difficult to establish a trend with ion- 
ization temperature. A more systematic approach will 
be presented in Pinto et al. (in preparation). 

15 www.sron. nl / files /HE A / SPEX / manuals / manual . pdf 



4.4. Comparison with atmosphere models 

The closest physical approach to derive global param- 
eters of the nova ejecta are non-LTE atmosphere models. 
Of particular interest for understanding nova evolution 
are the effective temperature and the mass of the under- 
lying WD. Here, we use the publicly availa ble tabulated 
non-LT E atmosphere models described by iRauch et al.l 
(|2010bD . The Tiibingen Model Atmosphere Package 
(TMAP) has been developed to study atmospheres of 
hot white dwarfs, where assumptions such as hydrostatic 
equilibrium and a plane parallel geometry are valid ap- 
proximations. In ^5.2.11 we discuss in more depth the 
limitation of these assumptions when applied to novae. 

We used all available data tables and obtained best fits 
with xspec (Arnaud 1996) for all four spectra extracted 
from the time intervals listed in Table [31 We carried 
out simultaneous fits to the separate RGSl and RGS2 
spectra and used the fluxed spectra only for plotting. We 
corrected for photoele ctric absorption us i ng th e TBnew 
module developed by (jWilms et al.|[2000l [20061 Wilms, 
Juett, Schulz, Nowak, in preparation) with a variable 
column density of neutral hydrogen, iVu, and abundance 
of neutral oxygen. 

The public database currently contains only models 
assuming log 5 = 9, covering a range of chemical abun- 
dances. In Fig. [Til the best fit models #007 and #003 are 
shown with blue solid and red dotted lines, respectively, 
in comparison to the spectrum from day 40.18. The top 
panel shows the entire spectrum, and the overall repro- 
duction of the data appears reasonable. The oxygen-rich 
Model #007 performs better at shorter wavelengths, sup- 
porting our earlier conclusion of a high O abundance. In 
the bottom panel, a zoom-in is shown, illustrating that 
some of the absorption lines are represented well. Values 
of reduced (number of degrees of freedom in brackets) 
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for model #007 are 3.7 (5511), 2.2 (5613), 21.8 (5616), 
and 4.5 (5504), for days 39.93, 40.03, 40.18, and 49.62, 
respectively. Note that the number of degrees of freedom 
is overestimated (and thus reduced underestimated) 
because a large number of parameters in the models are 
not counted as free parameters because they are not var- 
ied. Hard-coded values such as NLTE parameters have 
already been chosen in an early stage of model devel- 
opment and are not further considered during routine 
fitting, but they would likely change the results if mod- 
ified. Therefore, they are not necessa rily 'uninteresting' 
parameters in the sense of lAvnil (|1976f ) and may be of im- 
portance in deeper studies involving new computations 
for the special situation of an individual system such as 
V2491 Cyg. In particular, different values of logg = 9 
that are currently not publicly available have to be con- 
sidered, which will significantly increase the parameter 
space, likely yielding different results. 

The model parameters are listed in the right part of 
Table |3l where the ranges include the best-fit results 
from using all available models. The ranges do not re- 
flect the full accuracy of physical quantaties but only the 
ranges of parameters if different abundance classes are 
assumed. It is impossible to predict how the parame- 
ters would change if effects of expansion or extended ge- 
ometry were included. Radii were computed from the 
model normalization, R = 10~^^ x y^norm x d, with 
d = 10.5 kpc the distance, bolometric luminosities from 
the Stefan-Boltzmann law, and masses from the radii and 
logy — 9, M = R^g/G with the gravitational constant 
G= 6.674281 x 10"^ in cgs units. 

The parameter of highest interest is the effective tem- 
perature, for which we find the same trend with time as 
already found from the blackbody fits. With the result 
from the line profiles that the ejecta are extended, al- 
lowing us to view through different velocity layers, it is 
plausible that we are also viewing different temperature 
layers in the same spectra. Therefore, the definition of a 
single characteristic temperature does not have the same 
physical meaning as in a white dwarf evolutionary model 
and can not be used in the same way. We also caution 
that the underlying parameter log 17 = 9 is unrealistic for 
nova ejecta. The derived masses are completely unrealis- 
tic with some values much larger than the Chandrasekhar 
mass limit. We note that standard mass-radius relations 
do not apply for active novae since the radiative energy 
is not drawn from gravitational contraction but from nu- 
clear burning. Also, the bolometric luminosity computed 
from the Stefan-Boltzmann law, assuming the same ra- 
dius, is of order lO'^^ergs"^ and thus also unrealistically 
high for the late stage of evolution. 

It has to be kept in mind that the underlying model as- 
sumption is a plane parallel geometry, and strictly speak- 
ing, the calculation of a radius is in violation to this 
model assumption. Therefore, all calculations of WD 
masses and luminosities that have been obtained from 
hydrostatic plane parallel models have to be treated with 
greatest caution. 

With all these concerns in mind, these models repre- 
sent the data surprisingly well as can be seen in Fig. [11] 
making it tempting to trust this approach. Possibly, 
new data tables with lower log g may yield more realistic 
masses, however, spectral fits of models with more realis- 
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"Projected, =assumcd rest wavelength, assuming blue shifts of 
3300kms-i for V2491 Cyg, 2400 km s"! for V4743 Sgr, and 
1200kms-i for RS Oph. 



tic parameter combinations have not so far yielded better 
fits than can be seen in F ig. [TT] (see, e.g.. iNelson et al.l 
[2OO8I: IRauch"erani2010aD . 

We argue that a much more detailed discussion of re- 
finement of atmosphere models is required, which is be- 
yond the scope of this paper. While the data tables can 
be extended to include a larger range of logg, it will 
in our view be unavoidable that individual models are 
calculated for each nova, accounting for the individual 
conditions. Furthermore, more studies are needed to in- 
vestigate the effects on the derived effective temperature 
when accounting for the extended geometry and the ex- 
pansion. 

5. DISCUSSION 

5.1. Light Curves 

Based on the second half of the first observation, a pe- 
riodic signal of 2232.1 seconds (37.2 minutes, frequency 
0.448 mHz) may be present. We cannot entirely ex- 
clude red noise, but it is noteworthy that this signal 
is of the same order as the ~ 2500-s period detected 
bv lDrake etldl (I2003D in V1494 Aql or the 1300-sec pe- 
riod seen bv iNess et all (|2003D in V4743Sgr. It could 
be the s pin period of the WD or orig i nate from pul- 
sations ([Drake et all 120031: IKang et all I2006D . A de- 
tailed study of th e X-ra y light curves of V4743 Sgr by 
IDobrotka fc Nesd (|2010D revealed the presence of two 
nearby signals betwe en days 180 and 196 (see also 
iLeibowitz et"al|[2006l ). While roughly the same period 
was found to be persistent over the course of three years, 
small frequency changes were found by Dobrotka & NeH 

([Mo). 

During the SSS phase, the frequency changes scaled 
with changes in the blackbody temperature as pre- 
dicted by stellar pul sation theory (for details see 
IDobrotka fc Nesd l20ldl ). After day 526, a stable fre- 
quency was found in post-outburst X-ray light curves. 
However, in optical and UV light curves of V4 743 Sgr 
a slightly different frequency was found. They con- 
cluded that while pulsations may have been present dur- 
ing the SSS phase, the post-outburst X-ray period could 
be the WD spin period, modulated by pulsations dur- 
ing the SSS phase. Combined with the orbital period. 
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m V2491 Cyg (doy 40.18) 
H V4743Sqr (day 180) 
- Rest woveTength 

Expected wovelenqth 

if i>=-2400km/s 




28.5 29.0 29.5 30.0 30.5 31.0 

\ (A) 

Fig. 12. — Comparison of a calibrated Chandra/hFiTGS spec- 
trum of V4743 Sgr (day 180 after discovery; light blue) to the phase 
c spectrum of V2491 Cyg (starting day 40.18, red as in Fig.[3)l. The 
same lines as in Fig. [3 are marked, and the dotted lines indicate 
the expected blue shifts for an expansion velocity of — 2400kms~^ 
UNess et al.ll2003i) . 



V2491 Cyg (doy 40.18) m 
RSOpti (doy 39.7) X 0.40 m 
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if v=-1200l<m/s 
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Fig. 13.— Same as in Fig.[T2lfor RS Oph (light blue). Here, Une 
shifts of -1200kms-i are marked UNess et al.ll2007r) . 



iKang et alj (|2OO60 and lDobrotka fc NessI (|2010D conclude 
that V4743 Sgr is an intermediate polar. For V2491 Cyg, 
only two light curves are available, the first of which is 
interrupted by a dip that has to be excluded from the 
period analysis. In intermediate polars, the orbital pe- 
riod is usually significantly longer than the spin period 
of the WD, i.e. , there is no spin/orbit synchronization. 
iBaklanov et al.l ()2008() found a period of 2.3 hours for 
V2491Cyg, which is much longer than the 37.2-minute 
period, but it was not claimed to be the orbital period. 
The fact that a period can be seen in the X-ray light curve 
could imply that magnetic fields are involved, but not 
enough information is available to make a strong enough 
case for an intermediate polar on these grounds. Partic- 
ularly puzzling is that the period is not present in the 
second observation, while the X-ray light curve of an in- 
termediate polar should be modulated with the spin pe- 
riod of the WD. Spectroscopic evidence supporting the 
interpretati o n of an intermediate polar is discussed by 
iTakei et al.l (l20Tll) . 

5.2. Spectra 

5.2.1. Applicability of ttieoretical stellar atmospheres 

The ultimate goal of spectral analyses is to determine 
mass loss, composition, and kinematics. The require- 
ments for finding these parameters are well-exposed X- 
ray grating spectra that allow individual lines to be re- 
solved and sophisticated spectral models that include the 
requisite physics and reproduce the observed spectra. 

X-ray spectra of CCD resolution 100 eV) are 
not sufficient, as spectral models that are different in 
their details will be indistinguishable after folding them 
through the instrumental response of a CCD spectrome- 
ter. Too many different models can reproduce the same 
data, yielding non-unique results. Fine tuning of some 
parameters may improve the value of x^i but with the 
limited amount of information in the observed spectrum, 
the adjustments in the model are arbitrary. 

With the XMM-Newton and Chandra gratings, appro- 
priate spectra can be taken and are presented in this pa- 
per. The first requirement, well-exposed high-resolution 
spectra, has been met, however, spectral models need 
more development. The kinematics in the ejecta are rel- 
atively easy to determine from the blue shifts in the ab- 
sorption lines, but at the same time they pose the biggest 
challenge to self-consistent non-LTE atmosphere models. 

Another problem is that we depend on a very limited 
number of publicl y availabl e model s. An attempt has 
been made by R anch et al.l ()2010bf ) to make available 
an archive of TMAP models, covering a multi dimen- 
sional space of parameters. While the approach yields 
the most sophisticated models using a large database of 
transitions and treating radiation transport in full non- 
LTE, they are plane-parallel plus static and as such in- 
appropriate for nova ejecta while they are still expand- 
ing. We have shown in i)4.31 where Table |4] is discussed 
and Fig. [9] is shown as an example, that all photospheric 
absorption lines are blue shifted as well as significantly 
broader than the instrumental line broadening function. 
This can only be interpreted as expanding ejecta that are 
sufficiently extended to observe a range of layers o f dif- 
ferent velocity and consequently also temperature. iNessI 
((201Q) has shown that this phenomenon can also be ob- 
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served for other novae. The TMAP models account for 
non-LTE effects, however, the underlying assumptions of 
hydrostatic equilibrium and plane parallel geometry are 
far from reality for expanding nova ejecta, and our re- 
sults show that such simplifications are not justified. For 
more detailed discussion of the importance of the mis s- 
ing p h ysics, we refer to Hauschil dTet al.l ()1992l Il994bl lal. 
119951) ; iHauschildt fc Baron (^1993). Although, these ar- 
guments apply to an earlier phase of the nova evolution, 
the high-resolution spectra indicate that the arguments 
are the same. No models are publicly available that 
account for all the necess ary physics, and we refer to 
Ivan Rossum fc Nesi (|2010D who found that atmosphere 
models that account for the expansion yield vastly dif- 
ferent results from static models. 

Clearly, atmosphere models are required that can give 
a full account of the physics in the ejecta without neglect- 
ing the extended nature and/or expansion, including the 
effects of spherical and non-spherical geometry. Without 
such efforts, no accurate determination of a photospheric 
temperature is possible, even if the spectra can be well 
fitted. 

Meanwhile, blackbody fits are at most parameter- 
izations yielding no more than a spectral hardness 
parameter. Some limited, qualitative studies of spectral 
changes can be carried out, but no quantitative inter- 
pretation is possible. The advantage of blackbody fits is 
the small parameter space, yielding unique parameters, 
and blackbody temperature values respond sensitively 
to changes in spectral hardness. While different atmo- 
spheric processes may balance out to result in similar 
temperatures with the brightness decrease on day 40.03, 
it appears more likely that only the brightness has 
changed while the atmospheric structure remained the 
same. This would argue in favor of external reasons for 
the brightness variations. Meanwhile, the increase in 
blackbody temperature towards the second observation 
could be explained by either a higher photospheric 
temperature or changes in the atmosphere structure. If 
interpreted as an increase in photospheric temperature, 
this result would be consistent with a decreasing photo- 
spheric radius, exposing successively hotter material. 

5.2.2. Comparison to other novae 

V2491 Cyg is the third nova for which a bright SSS 
X-ray spectrum with atmospheric continuum and deep 
absorption lines was obtained in high spectral resolu- 
tion, after V4743 Sgr and RS Oph, and a comparison is 
instructive. 

V4743 Sgr 

In Fig. 1121 a comparison between a calibrated Chan- 
(ira/LETGS o bservation of the slower (^2 — 9 days, 
iMorgan et all [2003^ nova V4743 Sgr as observed on day 
180 after discovery (light blue) and the day 40.18 spec- 
trum of V2491 Cyg (red, see Table[3]) is shown. The spec- 
trum of V4743 Sgr is significantly softer than V2491 Cyg, 
indicating a lower photospheric temperature. While the 
lower E{B - V) = 0.25 for V4743S gr (as oppo sed to 
E{B -V) ^ 0.43 for V2491Cyg bv iRudv et al.ll200l 
could explain a softer X-ray spectrum, the Wien tail 
being at longer wavelengths makes a strong case for a 
lower effective temperature. In the bottom panels, the 
two wavelength ranges 22-27 A and 27-31 A are shown 



in detail. The same line labels as those used in Fig. [7] 
are included, but here, we assume a blue shift corre- 
sponding to — 2400kms~^, which is consistent with sev- 
eral photospheric lines in V4 743 Sgr such as Nvil at 
25 A ()Ness et al.|[200l . The line profiles in V4743Sgr 
are much more complex than in V2491 Cyg, some of 
them consisting of more than one velocity component. 
A similarly comp lex situation is encountered in KTEri 

ess et al.]l201Clt ). For example, the line profile around 
the Nvii line is extremely complex in V4 743 Sgr, while 
for V2491 Cyg, only two overlapping lines from N vii and 
Nvi are seen (see Fig. [9]). The interstellar Ni and Oi 
lines at 32.28 A and 23.5 A, respectively, are present in 
both novae at the same wavelengths. The O i line blends 
with the blue-shifted N VI7 line in V4743 Sgr. 

RS Oph 

More similar to V2491 Cyg is a grating spectrum taken 
39.7 days after the discovery of the 2006 outburst o f 
the recurrent, symbiotic nova RS Oph (|Ness et al.ll200l 
which evolved on a similar time scale (^2 = 7.9 days, 
pomisell et al. 2010). This is shown in Fig. [T51 where 
the calibrated C/iandra/LETGS spectrum of RS Oph (in 
flux units) was scaled by a factor 0.4. The observed 
blue shifts seen in the photosphe ric lines in RS Oph was 
-1200kms-i ijNess et all 120071 ) . and the dotted lines 
mark these projected velocities. 

The scaling factor of 0.4 is much larger than is required 
for th e different distances of 1.6 kpc for RS Oph (.Bodi 
IT987I ) and 10.5 kpc for V2491Cyg. Thus, at the respec- 
tive times of observation, V2491 Cyg was intrinsically a 
factor 17 brighter than RSOph. However, both novae 
were extremely variable around the respective times of 
observations. 

The spectrum of V2491 Cyg used for the comparison 
has been extracted from a phase of bright emission (day 
40.18) while the Chandra observation of R S Oph was 
taken during a time of comparatively low flux (i Ness et al.l 
I2007f l. In addition, there are uncertainties in the distance 
of V2491 Cyg, because the temporary secondary peak 
seen in V2491 Cyg (see bottom panel of Fig. [T|) disallows 
the use of the MMRD method. Under the assumption of 
the same luminosity at the respective times of observa- 
tion, the scaling factor between RS Oph and V2491 Cyg 
yields a distance of 2.5 kpc for V2491 Cyg. The soft 
tails of the X-ray spectra of V2491 Cyg and RS Oph are 
very similar, indicating the same amount of photoelec- 
tric absorption, which is consistent with t he measure- 
ment s of E{B -V) = 0.43 for V2491 Cyg (|Rudv et all 
120081 corresponding to iVn = 2.6 x 10^^ cm~^) and 
TVh = 2.4x 10^1 cm-2 for RS Oph (iHjellming et al.ll9 86). 

Noteworthy is a remarkable similarity in the continua 
between the X-ray grating spectra of V2491 Cyg and 
RS Oph (Fig. [13]) suggesting a surprisingly similar photo- 
spheric temperature at the times of the respective obser- 
vations. With CCD resolution, V2491 Cyg and RS Oph 
would be considered twins. However, at grating resul- 
tion, small but important differences can be seen, yield- 
ing information about the structure of the atmosphere. 
Most obvious are much smaller blue shifts and line widths 
of the photospheric lines in RS Oph, indicating slower 
expansion. Also, the Oi absorption edge at 23.5 A is 
slightly steeper in V2491 Cyg. While the Nvi and Nvii 
lines around 28.7 and 25 A reach about the same depth 
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in both cases, the O viii hne around 19 A is much deeper 
in V2491 Cyg. Since the same depth is reached for the 
hnes of two different ionization stages, Nvi and Nvii, 
the deeper Oviii hne cannot be a temperature effect, 
and a higher oxygen abundance in V2491 Cyg appears 
hkely. This could also explain the deeper O i absorp- 
tion edge. On the other hand, the Ovil hne at 21.6 A 
is only shghtly deeper in V2491 Cyg, but that could be 
explained by photoionization from O Vii to O VIII, which 
is supported by the presence of the O vii absorption edge 
at 16.7 A in V2491 Cyg (see Fig.® that seems less pro- 
nounced in RS Oph. 

Unidentified lines 
While the optical evolution was different, yielding 
a rebrightening in V2491 Cyg that was not seen in 
RS Oph, V2491 Cyg wa s suggested t o be a r ecurrent nova 
(IPagnotta et al.l l2009t IPage et all 120101 IRagan et all 
I2010D . and the spectral similarity to one of the most 
famous recurrent novae, RS Oph, may be an additional 
indicator. 

While the measurement of line shifts relies on well- 
identified lines, a number of absorption lines cannot be 
identified. The comparison with other novae gives us 
an idea of the formation conditions of unidentified lines, 
aiding future identification efforts. At 25.85 A, a broad 
line is seen in V4743Sgr that coincides with a line at 
26.06 A that is marked as unidentified in Fig.[71 With the 
lower photospheric temperature of V4743 Sgr, this line 
could be formed at relatively lower temperature since it 
is much more pronounced in V4743 Sgr. Another uniden- 
tified line in Fig. [71 with a projected rest wavelength of 
29.45 A, is also present in V4743Sgr with roughly the 
same width and strength. The only difference is the line 
shift, which is similar to the Nvi line at 28.8 A. This line 
might be formed at a slightly higher temperature than 
Nvi which is deeper in V4 743 Sgr. The line at 30 A pro- 
jected rest wavelength is not present in V4 743 Sgr and 
could be formed at temperatures that are higher than 
in V4 743 Sgr. The unidentified lines are listed in Ta- 
ble [5] with projected rest wavelengths in the first column 
and observed wavelengths for V2491 Cyg, V4743 Sgr, and 
RS Oph (see below) , respectively. 

The unidentified lines with projected rest wavelengths 
26.06, 29.45, and 30.0 A are also present in RS Oph. The 
26.06-A line is weaker in RS Oph than in V4743Sgr, but 
stronger than in V2491 Cyg, while the 29.45-A line has 
about the same strength. This supports the conclusion 
that the 26.06-A line is a low-temperature line while the 
29.45-A line is formed at higher temperatures. In con- 
trast to V4743Sgr, the 30.0-A line is clearly present in 
RS Oph and is deeper than in V2491 Cyg. Additional 
unidentified lines are seen in RSOph, e.g., at 25.6, 25.7, 
26.1, 27.3, 27.5, 27.6, 28.0, and 30.3 A, but these lines 
are not seen in V2491 Cyg and V4743 Sgr. Two of the 
unidentified lines are seen in all three novae at wave- 
lengths yielding the same blue shifts as other photo- 
spheric lines. Meanwhile, a large number of lines are 
only present in the spectrum of RS Oph, arising in the 
same wavelength range as a pre-SSS emis s ion line com- 
ponent around day 26 ()Ness et al.l l2009bl : iNelson et al.l 
[2008.1 . 

The emission lines in the RS Oph spectrum at 28.8, 



29.1, and 29.5 A originate from the Nvi He-like triplet. 
These and other components of the red wings of absorp- 
tion lines such as Oviii (18.97 A) and Nvii (24.8 A) 
are related to th e P Cygni-like profiles discussed by 
iNess et al.l (|2007D . They could partially originate in 
shocked plasma and might not be related to the atmo- 
spheric emission. 
Other novae 

Other nova SSS spec tra such as those of V1494Aql 
(jRohrbach et al.l 120091) . are more complicated and/or 
less well exposed. AU high-quality SSS spectra show 
significant blue shifts in their absorption lines, and the 
conclusion is un avoidable t hat expansion continues into 
the SSS phase (|Nes£ll2010D . It is noteworthy that the 
expansion velocities can be quite different, yielding the 
lowest expansion velocity in RS Oph, even though the 
spectral shape of the continuum and the depths of most 
absorption lines are remarkably similar to V2491 Cyg. 
The lower expansion velocity in RS Oph may be ex- 
plained by the presence of a dense stellar wind of the 
companion that has absorbed some of the kinetic energy 
during the earlier phases of the outburst. This results 
in significant deceleration occur ring by ~ 40 days after 
outburst (e.g. iBode et al]|2006[ ). While some evidence 
for the compani on in V2491 Cyg to be an evolved object 
is presented by IRibeiro et al.l ([201 II ) , the mass loss rate 
is likely to be much lower than in RS Oph. Meanwhile, 
in V4 743 Sgr, some of the absorption lines contain at 
least three velocity components while V2491 Cyg and 
RS Oph contain only one velocity component. Since 
V4 743 Sgr is a slower nova, it may be possible that the 
ejecta in slower novae are generally more structured than 
those of faster novae. This would be con sistent with 
modehng of comm on envelope interaction ([Llovd et al.l 
119971 : [Porter et al.|[l9981 , but other slower novae need to 
be studied for firmer conclusions. 

5.2.3. High-ampUtude vanattons 

During the earlier of the two XMM-Newton observa- 
tions, lasting from day 39.93 to 40.37, V2491 Cyg was 
highly variable in X-rays, including a three-hour dip. 
Several novae have been observed with an episode of 
high-amplitude variations during their early SSS phase, 
but an explanation is still pending. The first observed oc- 
currence of a sudden decay in the X-ray output of a nova 
was encountered in a Chandra ob servation of V4743 Sgr 
(see Fig. [H and INess et~alll2003f ). The count rate de- 
clined almost to zero, leaving behind an emission line 
spectrum, originating from photoionization plasma with 
some emission lines that were rapidly declining (Ness et 
al. in prep). At that time, no monitoring was possible, 
and the nova did not rebrighten during the same obser- 
vation. Two we eks later, the nova was bright again, but 
highly variable ([Leibowitz et al.ir2006[) . The dip in the 
first XMM-Newton observation of V2491 Cyg occurred 
on the same time scale as the decay in V4743 Sgr, as il- 
lustrated in Fig. 1141 and the same phenomenon may have 
been observed. 

Dense Swift monitoring of RS Oph revealed multiple 
declines ([Osborne et al.l [201111 which occurred on time 
scales that are also similar to V2491Cyg (Fig. . This 
phenomenon has not been fully explained yet. One 
suggestion is temporary photospheric expansion, which 
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Fig. 14. — Comparison of the X-ray light curve of V2491 Cyg after 
day 39.6 with rescaled Chandra and Swift light curves of V4743 Sgr 
with the de cay on day 180.6 and of RS Oph with one of the minima 
reported bv lUsborne et all l |2011 ). respectively. 

would be supported by the disappearance of periodic os- 
cillations during the dip which likely originate from close 
to the WD. Such deep dips could also mean that the 
concept of constant bolometr ic lum inosity is wrong as 
suggested by Starrficld ct al] ()2009( ). However, intrinsic 
changes in bolometric luminosity are not necessarily eas- 
ier to understand than changes in the structure of the 
ejecta. A change in luminosity would imply a change 
in the ionization parameter, ^ = L/{nr'^), of the photo- 
ionized gas, with a subsequent evolution of the photo- 
ionization balance. Since the density of the photo- ionized 
absorbing gas, n, and distance to ionizing source, R, are 
unlikely to change significantly, intrinsic changes of the 
bolometric luminosity would have lead to equally sig- 
nificant changes in the ionization parameter. A detailed 
study using photo-ionization models is in progress (Pinto 
et al. in prep). 

6. SUMMARY AND CONCLUSIONS 

A small sample of extremely well exposed grating spec- 
tra of novae during their SSS phase has been gathered. 
The fast classical nova V2491 Cyg was observed twice 
during its extremely short SSS phase, covering two dif- 
ferent phases of evolution. The first observation was 
taken early while the X-ray brightness was highly vari- 
able, containing a deep dip, while the second observation 
was taken during the decline with a lower count rate and 
less variability. The dip evolved on a ti me scale similar to 
the decay in the slower nova V4743 Sgr (INess et al.ll200l 
and the variations du ring the early SSS phase in RS Oph 
(jOsborne et al.|[20Tl|) . High-amplitude variations seem 
to be a common phenomenon that has now been observed 
in many novae. 

The contemporaneous UV light curves show short-term 
variability, but no UV data were taken during the dip in 
the first X-ray light curve. A 37.2-minute period was de- 
tected in X-rays before and after but not during the dip. 
It could reflect the WD spin period, but the statistical 
significance is not high. 

The large majority of absorption lines in the high- 
resolution RGS X-ray spectra yield a single velocity 
component with blue shifts ranging between —3000 and 
— 3400kms~^. Additional interstellar lines are found at 
their rest wavelengths. A few deep absorption lines are 
currently not identified. The same lines are also present 
in the SSS spectra of V4743 Sgr and RS Oph at a differ- 
ent blue shift, indicating that they are of photospheric 
origin. 

The narrow range of blue shifts is different from 
V4743 Sgr where at least three velocity components can 
be seen. This could be explained by the slower evolution 



of V4743 Sgr. The SSS spectrum of RS Oph is remark- 
ably similar in the shape of the continuum, yet is very 
different in the properties of the absorption lines, yielding 
much smaller blue shifts. The lower velocities in RS Oph 
may be due to the presence of the pre-outburst red giant 
wind, while the higher blue shifts in V2491 Cyg could be 
the explanation for the short duration of the SSS phase. 

Archived atmosphere models yield good fits to the 
SSS spectra, but the model parameters are unphysical. 
This can be explained by the simplifications in the 
model, assuming hydrostatic equilibrium and a compact 
atmosphere without expansion. 

The conclusions are the following: 

1. Significant blue shifts in all photospheric absorp- 
tion lines demonstrate that the nova ejecta are ex- 
panding. 

2. Significant line broadening of most photospheric 
absorption lines demonstrates that the ejecta are 
extended, exposing a range of velocity layers and 
thus also temperature layers. 

3. Our conclusions from line profile analysis indicate 
that publicly available stellar atmosphere models 
do not contain enough physics for nova ejecta, even 
though they yield surprisingly good fits to the data. 

4. Comparison with other novae of lower and higher 
photospheric temperatures allows us to spot so far 
unknown nova absorption lines and associate them 
to higher or lower temperature environments. 

5. Remarkable spectral similarities between 
V2491Cyg and RS Oph indicate similar system 
characteristics. Assuming the same luminosities, 
the distance to V2491Cyg would only be 2.5 kpc. 

6. The different properties of the absorption lines be- 
tween V2491 Cyg and RS Oph indicate that similar 
global parameters can still lead to a different at- 
mospheric structure. 

7. Slower novae such as V4 743 Sgr may contain more 
complex absorption line systems than faster novae 
such as RSOph or V2491Cyg. 

8. The single dip in the X-ray light curve occurred on 
the same time scale as in V4743 Sgr and RS Oph, 
and the underlying process may therefore be com- 
mon to all novae in which an early variation phase 
is observed in X-rays. 
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